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Consider a system consisting of n d-dimensional quantum 
particles and an arbitrary pure state |$) of the whole system. 
Suppose we simultaneously perform complete von Neumann 
measurements on each particle. One can ask: what is the 
minimal possible value Sfif] of the entropy of outcomes joint 
probability distribution? We show that S[^f] coincides with 
entanglement entropy for bipartite states. We compute S[9] 
for two sample multipartite states: the hexacode state |H) 
(n = 6, d = 2) and determinant states [Det n ) (d = n). The 
result is S[H] = 4 log 2 and S"[Det n ] = log(n!). The general- 
ization of determinant states to the case d < n is considered. 

PACS Nos. 03.67.-a, 03.65.Bz 



I. INTRODUCTION AND MAIN RESULTS 

Quantum information theory has many interesting fea- 
tures which have no classical analogue. One of them is 
entanglement or quantum correlations. It has been the 
object of intensive study for last years because it is the 
entanglement that makes possible to develop effective al- 
gorithms solving many tasks in computing, communica- 
tion, and cryptography, see for example Q and references 
therein. 

However at the present moment the canonical defini- 
tion of entanglement is missing and the question how to 
quantify the degree of entanglement in a given multipar- 
tite quantum state remains open. The only reasonable 
constraint on a functional which pretends to be an entan- 
glement measure is the monotonicity under certain class 
of local quantum operations Q . 

An important example of such functional is entangle- 
ment entropy ||. For a pure state \^) of bipartite system, 
its entanglement entropy E[^] is defined as 



The functional E^] has very simple physical sense. 
Suppose that each of two parties, between which the state 
1 5') is distributed, performs complete von Neumann mea- 
surement on his part of the system. Such joint measure- 
ment is a complete measurement on the whole system. 
Its outcome is a random variable whose probability dis- 
tribution depends upon the pair of complete von Neu- 
mann measurements chosen by each of two parties. This 
choice is equivalent to the choice of orthonormal basis in 
each party Hilbert space of states. One can ask: what 
bases should be chosen by each party to minimize the en- 
tropy of the outcomes joint probability distribution and 
what is the minimal value of this entropy? If the state 
|*) is factorizable, i.e. |W) = (g) |*2), then the an- 
swer is trivial: the i-th party should complement | ) 
to a complete basis by any way and this choice yields 
zero entropy because the measurement outcome will be 
"(^1,^2)" with the probability one. If 1*) is entangled, 
one can easily show (see Section ||) that the z-th party 
should perform measurement in the basis where its den- 
sity matrix pi is diagonal (the Schmidt basis) and the 
minimal entropy of the outcomes coincides with £ I [\I']. 
Thus we can interpret entanglement entropy E[^] as the 
minimum of outcomes entropy over all choices of local 
complete von Neumann measurements. 

If we consider entanglement entropy from this point of 
view, it can be naturally defined for arbitrary multipar- 
tite states. Suppose a system consists of n d-dimensional 
quantum particles distributed between n remote parties 
and let |\I>) S (C")® n be arbitrary pure state of the whole 
system. Denote Bi the orthonormal basis in C d chosen by 
i-th party and let \Bi(j)) £ C d be the j-th basis vector in 
the basis Bi, where i € j S [l,d\. Orthonormality 

condition implies that (Bi(j)\Bi(j')) — 6jj>. Now let us 
define the functional S[$f] according to: 



E [*] = ~ £ Pi kga ft 



(1) 



S[9] 



inf S[9,B U 

Bx,...,B„ 



,B n 



where pj's denote the eigenvalues of reduced density ma- 
trices pi = tr 2 (|\I>)(*|) and p 2 = tri(|*)(*|) (they have 
the same spectrum). This particular measure is distin- 
guished because in the asymptotic limit (i.e. when one 
takes a large number of copies of a given shared state) 
any monotonic functional of bipartite state up to triv- 
ial rescaling coincides with entanglement entropy, see || . 
Unfortunately, when the system is divided to three or 
more local parts, entanglement entropy is not defined. 



S[^,B 1 ,...,B n } = -J2p(j)^g 2 p(j), 



P (j) = p( 3l , ...,j n )= K*|Bl0l), ■ ■ ■ , Bn (jn)) I 2 , 



(2) 



(the sum over multi-index j is the sum over all possible 
iii • • • i3n G [lj d]). It tells us to what extent the parties 
may decrease the entropy of the outcomes distribution 
by varying the bases in which they perform the measure- 
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ments. As was said above, S[^f] = E[^>] for bipartite 
states, so we will call S[*ff] as entanglement entropy. Its 
properties immediately following from the definition are 
< < nlog 2 d; S[^f] = iff the state |#) is factor- 
izable; ® <&] = S[^] + (here we mean that |\I>) 
and |$) are the states shared by n and m parties, while 
|\& (g> 4>) is shared by n + m parties); S[$f] is continuous 
functional of the state |^). 

Computation of entanglement entropy in the multipar- 
tite case n > 2 is a difficult task. Probably, for generic 
quantum state it can be solved only numerically (note 
that the number of parameters to be optimized in the 
definition (|J) grows as 0{nd 2 )). It is relatively easy to 
get the upper bound on S\9] — one just needs to choose 
tentatively some basis for each party. As the lower bound 
on S\*S>] one can take von Neumann entropy of the mixed 
state of any group of parties (see Section [II). However 
for generic state this lower bound is too weak. It can be 
improved if the state has some symmetry. 

In this paper we consider two types of symmetry on 
example of determinant states |Det„) £ (C")® n , see [[|, 
and six qubit hexacode state |H), see ||. The determinant 
state |Det„) is invariant under unitary transformations 
U® n where U £ SU(n) is arbitrary one-party unitary op- 
erator. The hexacode state is in some sense 'maximally 
uniform' pure state — if we divide six qubits into two 
equal groups by arbitrary way then the mixed state of 
each group will be absolutely uniform. Due to these spe- 
cial properties, entanglement entropy of the states |Det„) 
and |H) can be exactly computed. 

Another reason for our interest to these particular 
states is that their entanglement entropy is rather close to 
the upper bound n log 2 d, so they are near-maximally en- 
tangled states. For the hexacode state (n = 6, d = 2) the 
computation yields S[H] = 4, see Section VI, while for 
determinant state ( 



d) one gets SpeW] = log 2 (n! 



see Section [V. For large n we can approximately write 
SfDetn] m [1 — 1/ ln(ri)] n log 2 n. It means that deter- 
minant states asymptotically saturate the upper bound 
for normalized entropy: linin^oo S[Det n ]/n log 2 n = 1. 
In Section [v| we construct the generalized determinant 
states |Det„ ;d ) £ (Q d )® n defined if n = pd p , where p 
is arbitrary integer, such that SpDet,^] = log 2 [(eF)!]. 
For fixed d and large n (i.e. large p) we can write: 
S'petn^] ~ [1 — 1/ ln(n)] n log 2 d which again saturates 
the upper bound for normalized entropy. Note that the 
factor [1 — l/ln(n)] grows sufficiently slow, e.g. it is 
equal 0.9 for n = e 10 » 2 • 10 4 . We will see that 
for qubits (d — 2) the tensor powers of hexacode state 
|H® m ) have entanglement entropy greater than determi- 
nant states |Det„ 2) h" the number of qubits n < 60. 

We make concluding remarks and discuss some open 
questions concerning the generalized entanglement en- 
tropy in Section VII . Most interesting question concerns 
the stability of the definition (|^) under the extension of 
each party space of states. 



II. BIPARTITE SYSTEM 

It is known that up to local unitary operators, any 
state |W) € C d <£> C d of bipartite system is specified by 

its Schmidt coefficients {pi}i=i...d, Pi > 0; Hi=iPi = !• 
Being invariant under local unitaries S[*£] is a functional 
of the Schmidt coefficients only. Thus one suffices to 
compute S[&] only for the special states 



I*) 



'Pi m; 



(3) 



where {\i} £ C d }i = i...d is the standard basis of C d . If in 
the definition (Q) we tentatively choose the standard basis 
for both parties, i.e. \B\(i)) — \B 2 {i)) — \i), i £ [i-,d], 

then S[$>, 61,62] = -J2t=iPi^°S2Pi an d thus we get 
< £?[*], see (0). We can also prove that E[&] is 
simultaneously the lower bound for S[^f]. Indeed, let 6*, 
62 be the optimal choice of the bases (i.e. such that 
5[*] = S[^, 61,62]). Consider the density matrix p x of 
the first party only: p x = tr a (|*)<*|) = E*=iJ*l*)<*1- 

Denote Pl (i) = (BJ(0IPi|B*(i)> = E^K^WIi)! 2 
the distribution of the first party outcomes in the optimal 
basis. Because the entropy of partial distribution can not 
exceed the entropy of joint distribution, we have: 



S[*]>-5>i(*)log 2 [pi«]- 



(4) 



i=l 



Using the concavity of the function —x log 2 x and normal- 
ization 2j=i K^i (01?) 1 2 = 1' we § e ^ the next estimate: 

d 

-pi(i)log 2Pl (i) > - ^|(6 1 *( l )|j)| 2 felog 2 p J ). (5) 
j=i 

The summation over i can be carried out taking into 
account the normalization Ei=i l(^i(Ulj)| 2 = !• Thus 
from we can infer S[$>] > E[ty] and consequently 



(6) 



III. 



CONNECTION WITH VON NEUMANN 
ENTROPY 



Suppose the system consists of n <i-dimensional quan- 
tum particles distributed between n remote parties and 
let |W) £ (C d )® n be arbitrary pure state of the whole sys- 
tem. Let us choose a group X of k parties, for example 
X = {1,2, . . . ,k}. The chosen group of parties shares the 
mixed state p x = trj^ x (\^)(^\)- Suppose 6* is the opti- 
mal basis for the i-th party. Denote p x ... ,ik) = p x (i) 
the optimal outcomes distribution for the parties from X 
only: 
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p x (i) = (Bl(h), . . .^KiMBKh), . . .,B* k (i k )). (7) 

Let S\p x (i)] be the entropy of distribution (R). By repeat- 
ing the arguments presented in Section Owe can show 
that 



tr(p x log 2 p x ) < S[ Px (i)} < 



(8) 



Thus von Neumann entropy of the mixed state p x can 
serve as the lower bound on entanglement entropy. Of 
course the group of parties X can be chosen by arbitrary 
way. 

Note that the density matrix of the parties which were 
not selected to X has the same (positive) spectrum as 
p x . It means that one suffices to consider the groups 
of A; < n/2 parties and the best lower estimate on S[$>] 
which we can hope to achieve is (n/2) log 2 d. 

As a good example, consider three qubit GHZ state 
| GHZ) = 2-i(|0,0,0) + |1,1,1)), d = 2, n = 3. The 
density matrix of the first qubit is absolutely uniform: 
pi = (1/2)1 Thus 5[GHZ] > -tr(pi log 2 px) = 1. On 
the other hand we can choose tentative bases B\ = B2 
= B3 = {|0), |1)} which provide us with upper estimate 
5[GHZ] < S[GRZ,B 1 ,B2,B 2 ] = 1, so that 5[GHZ] = 1. 



IV. DETERMINANT STATE 

Let us consider the multipartite system with d = n. 
Choose the standard basis {\i) S C n }i = i... n in each copy 
of C" and consider the state |Det„) G (C n )®" defined as 



|Det„) = (n\y 



(9) 



where ^ ^ is completely antisymmetric tensor of the 
rank n and the sum is over all «i, . . . , i n £ [1, n]. This 
state was used in Ref. || to study the limitations on the 
pairwise entanglement in multipartite systems. We will 
call the family {|Det„)} n the determinant states. Note 
that |Det 2 ) is EPR singlet state 2 - 4(|l,2) - |2,1)). The 
purpose of this section is to prove the formula 



S[Det„] =log a (n!). 



(10) 



This result is immediate consequence of the following 
property of determinant states: 



sup 

pi\<ln)=l, ie[l,n]. 



|(Det r 



(n!)" 1 , (11) 



(we employ standard designation \<j>i, . . . , <f> n ) = |</>i) 
® ■ ■ ■ ® |^n))- It tells us that the projection of |Det„) 
on any factorizable state has the norm at most (n!) . 
To prove (|ll|), we first note that the state |Det„) is SU(n) 
singlet, i.e. for any V £ SU(n) we have ^^"iDetn) 



= |Det„). Therefore, if the projection of |Det„) on the 
state \(f>i, . . . , 4> n ) is the highest one, then the projection 
of |Det„) on the state \V<f>i, ■ ■ ■ , V<j) n ) is also the highest 
one. So while looking for the maximum in ( |ll| ) we can 



fix one of 



e.g. put 



But according to 



definition (|) we have: 
|(Det n |0i, . ..,</> n _i, n) | 5 



- 1 (Det n _i \<t> x ,. 



(12) 



the em- 

l\®(n-l) 



Here by abuse of notations we denote |Det n _i 
bedding of determinant state |Det„^i) £ (C n 
into the space (C n )®(™ -1 ) (the space C" -1 is embedded 
into C n by adding zero n-th component to all vectors). 
Although in (|l|) \(f>i) £ C n , the righthand side of (|l|) 
achieves the maximum when all states |0i), . . . , |</>„_i) 
have zero n-th component. It implies that 



1 



sup 

(^|0 4 )=l, i£[l,n]. 

sup 

\<Pi)ec n -\ 

(0,10,) = !, ie[l,n-l]. 



|(Det„|0i 



|(Det„_ 



which by induction leads to equality ( [Tl|) 



Now let us explain why ( |TT| ) implies (Q) . Suppose B* 
is the optimal basis for the z-th copy of C n , i £ [l,n], 
i.e. S[Det n ] = S[T)et n ,Bt,...,B*}. Let p*(i u . . . , i n ) 
= \(Bl(ii), . . . ,£>*(i„)|Det„)| 2 be the optimal distribu- 
tion. According to (pT|), p*(ii, . . . ,i„) < (n!) _1 for any 
outcomes i±,...,i n and thus S'petn] > log 2 (n!). On 
the other hand, we can tentatively suggest all parties 
to perform the measurements in the standard basis, i.e. 
\Bi(j)) = i,3 G [!,»]■ Then 5[Det n , Bi, . . . , B n ] 
= log 2 (n!) which tells us that S^Det^ < log 2 (n!) and 
thus that S^Detn] = log 2 (n!). 



V. GENERALIZED DETERMINANT STATE 

Suppose now that the dimension d of each particle is 
fixed. If the number of particles is 11 = pd p for some in- 
teger p, the space (C d )® n can be identified with (C d ")® dP 
and thus determinant state (Det^p) has its counterpart 
in (C rf )®". This simple observation allows to construct 
the state |Det„, d ) £ (C rf )®", n = pd p such that 



S[Det r , 



log 2 [(rf 



pdP. 



(13) 



Note that although |Detdp) and |Det nj( j) represent one 
and the same state, SfDetn^] might be greater than 
S^Detdp] because in the first case we have less freedom in 
the choice of bases in (||). 
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Let us explain the construction of the state |Det„ i( j} 
and derive ( |l3| ) on example of the qubits, i.e. d = 2. 
Consider any one-to-one map (p which maps the integers 
on the interval [1,2 P ] to binary strings of the length p, 
e.g. 



p(2) 
<p(2 p - 1) 



(0,0,. 
(0,0,. 

(1,1,. 
(1,1, 



,0,0), 
,0,1), 



1.0) , 

1.1) . 



(14) 



Define the state |Det„, 2 ) G (C 2 )®", n 
~ ^2 e ii,-,hp\ ( P(h),- 



Det, 



: p2*> as 



(15) 



where the sum is over all ii,...,i 2P G [1,2 P ] and we 
omit the normalizing factor [(2 P )!]~2. It is written here 
in the standard qubit basis {|0), |1)}. While computing 
SfDet^a] we minimize over the choice of p2 p one-qubit 
bases, see (||). But any such choice is also the choice of 
2 P bases in 2 P copies ofC 2 ". Therefore we can say that 



5[Det n , 2 ] > S[Det 2P ] = log 2 [(2")!]. 



(16) 



On the other hand, if we will tentatively measure each of 
p2 p qubits in the basis {|0), |1)}, the entropy of the out- 
comes distribution will be exactly log 2 [(2 p )!], see (15). 
Therefore £?[Det„ )2 ] < log 2 [(2 p )!] and thus the equal- 
ity ( |l3|) is proven (the proof for arbitrary d copies the 
proof for d — 2). 

As was already mentioned in Section the determinant 
states are near-maximally entangled states if the num- 
ber of parties n is sufficiently large. As an illustration 
let us consider determinant states |Det„ j2 ) correspond- 
ing to n qubit system. In Table I we present the normal- 
ized entanglement entropy S^Detn^J/n and the number 
of qubits n for p = 1, . . . , 5, and p = 10. 

TABLE I. Normalized entanglement entropy of determi- 
nant states |Det n ,2)- Here n = p2 p is the number of qubits. 



p 


n 


5[Det„,2]/n 


1 


2 


0.50 


2 


8 


0.57 


3 


24 


0.64 


4 


64 


0.69 


5 


160 


0.74 


10 


10240 


0.86 



Table I suggests that we could try to find some state of 
O(l) qubits which has entanglement entropy greater than 
the determinant states. This is the purpose of the next 
section. 



VI. HEXACODE STATE 

The hexacode state was originally defined in the con- 
text of quantum error correcting codes. It was associated 
with certain maximal self-dual linear subspace of GF(2) 6 , 
see ||, p. 30 for details. In this Section we present the 
alternative and more explicit definition of this state only 
briefly discussing its connection with quantum codes. We 
also prove the equality 

5[H] = 4, (17) 

announced in Section ffl. 

Consider graph G = (V, E) shown on Fig. 1 with the 
set of vertices V = {1,2,3,4,5,6} and the set of edges 
£={(12), (13), (14),... ,(56)}. 




FIG. 1. Graph G used in the definition of the hexacode 
state. 

We associate a qubit with each vertex i G V. Let Aij be 
6x6 adjacency matrix of G, i.e. Ay = 1 if (ij) S E and 
Aij = if (ij) E, — Aji. The diagonal elements 
An will not appear anywhere. Then six qubit hexacode 
state |H) G {C 2 ) m is defined as follows: 



|H) = (2 6 )-* J2 (-l) a{x) \x), 

xGB e 

o,{x) — ^ ' AjjXjXj , 



(18) 



where B 6 denotes the set of all binary strings of the 
length 6 and \x) = \x\, . . . ,x§). 

Note that |H) can also be defined in terms of stabilizers 
operators. Denote af the Pauli matrix a a acting on the 
«-th qubit and assign to each vertex of the graph the 
operator 



Xi=ot n 



o* v iev. 



(19) 



They commute with each other and stabilize the state 
|H), i.e. Aj|H) = |H). Six operators Xi generate the 
group S of all stabilizers, |5| = 2 6 . Each stabilizer from 
S is a tensor product of several Pauli matrices (proba- 
bly with '-' sign). As was shown in Q, any nontrivial 
stabilizer from S is tensor product of at least four Pauli 
matrices, so that |H) is an additive quantum code cod- 
ing qubits into 6 qubits with the minimal stabilizer 
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weight 4. Of course, such quantum code can not be used 
to protect quantum information from the errors. It is 
just symplectic state with some special properties. Note 
however that if some symplectic state has the minimal 
stabilizer weight d and if any [(d — l)/2] or less qubits 
were decohered then the syndrom measurement allows 
to determine the positions of the decohered qubits. 

The proof of the equality (|l7|) consists of three steps. 
The first step is to prove that S[H] < 4. On the second 
step we establish remarkable symmetry of the state |H) 
which is used on the third step to prove that S[H] > 4. 

1) Let us tentatively choose the bases Bi = {|0), |1)} 
for the qubits i = 2,3,4,5 and B x = B 6 = {|+),|-)} 
where | ±) = 2~2 (|0)±|1)). Simple calculations show that 
the distribution of the outcomes p{x) = p(xi, . . . , xq) 
= \(Bi(xi), . . . ,Be{xe)\H)\ 2 measured in these bases is 
following: 



p(x) 



j_ ^ i ii=i 2 + i 3 +i 4 (mod2), 
1 xe — xs + X4 + x^ (mod 2), 



0. otherwise. 



(20) 



Thus 5[H, Bi, . . . ,B 6 ] = 4 and consequently 5[H] < 4. 

2) Suppose the vertices of G are colored by black and 
white colors such that there are three black and three 
white vertices. Denote B and W the subsets of black 
and white vertices, B U W = V. Consider three-qubit 
density matrix p w describing the state of white qubits 
only: p w = trs(|H)(H|). We claim that for any partition 
V = BU W, p w = 2~ 3 1 where 1 is the unital matrix. In 
other words, the mixed state of any triple of the qubits 
is absolutely uniform. To verify this property, fix some 
partition and renumber the vertices of the graph to make 
W = {1,2,3} and B — {4,5,6}. The adjacency matrix 
A then can be split to four 3x3 blocks: 



A 



(21) 



which are adjacency matrices between white and white, 
white and black, black and white, black and black ver- 
tices (A wb = (A bw ) T ). The density matrix p w depends 
only upon A ww and A bw . Simple calculations yield: 



. 2 - 6 £ ( _ ir ^ 

2GB 3 



y 

'(*+y) 



(22) 



where x = (xi,x 2> x 3 ), y = (y 1 , 3/2,2/3), z = (z 1 ,z 2 ,z 3 ) 
and we treat A ww , A bw as 3x3 matrices over binary field 
acting on binary vectors. Observe that the sum over z 
is zero if A bw (x + y) ^ (0, 0, 0) and is equal to 2 3 if 
A bw (x + y) = (0,0,0). One can explicitly verify that 
the matrix A bw is nondegenerate over binary field M for 
any partition V = B U W (note that due to the sym- 
metry of the graph G, there are only three nonequiv- 
alent partitions, so this verification is very simple). It 



y. There- 
= 2~ 3 , i.e. 



means that A bw (x + y) = (0, 0, 0) only if x 
fore (x\p w \y) — if x ^ y and (x\p w \x) 
p w = 2~ 3 1, see also ||. 

3) Suppose the optimal basis for the i-th qubit is 
B*, i.e. S[B] = S[R,B*,...,B* 6 ]. Let p*(x u . . . ,x 6 ) 
= \(B*(xx), . . . , Bg(x 6 )\R)\ 2 be the optimal probability 
distribution. Consider any partition V = B U W. The 
probability distribution of three outcomes {xi : i G W} 
measured at white vertices only is 



^2 P*(x 1 ,...,x 6 ) 



Xi : i£B 



W B*(Xi)\p w \ <S>idW B*(Xi)) = -, 



(23) 



regardless of configuration {xj : i £ W}. In other words, 
the optimal distribution has a nice property: the par- 
tial distribution of any triple of bits is absolutely uni- 
form, see ||. Call such property of the distribution as 
3- uniformity. An example of 3- uniform distribution is ab- 
solutely uniform distribution. There are also 3-uniform 
distributions which are not absolutely uniform, e.g. the 
distribution ( po|) . Denote Pg the set of all 3-uniform dis- 
tributions of six bits. In Appendix ^| we show that 



inf S[p(x)} 
v(x)ePg 



(24) 



where S'[p(x)] is the entropy of probability distribution 
p(x). We know that p*{x\, . . . , Xq) £ P| and thus we 
have 5[H] = S\p*(xi, . . . ,x§)] > 4. This completes the 
proof of (0 



The definition like (18) can be used to assign a state 
|G> S (C 2 )®l y l to any unoriented graph G = {V, E). Rea- 
soning as above, one can show that if |V| = 2m and the 
adjacency matrix A bw is nondegenerate over binary field 
for any partition V = B U W, \B\ = \W\ = m, then 
in the state |G) any m qubits have absolutely uniform 
density matrix 2 _m 1. Such state |G) can be called max- 
imally uniform pure state, because any subset of qubits 
which is not forbidden by Schmidt constraint to have ab- 
solutely uniform density matrix do have absolutely uni- 
form density matrix. This property of the quantum state 
is interesting by itself. Surprisingly, for to < 15 (i.e. 
for |V| < 30) appropriate graphs exist only for to = 3 
(e.g. the graph shown on Fig. 1) and for m = 1 (e.g. 
V = {1, 2} and E = {(12)}). It follows from the bounds 
on additive quantum codes. Indeed, consider a state 
|G) assigned to such graph. Like as hexacode state, we 
can specify |G) by stabilizer operators (09) which gen- 
erate the group S of stabilizers of order 2^ m . Any sta- 
bilizer X £ S is a tensor product of several Pauli ma- 
trices (possibly with a sign '-') and X|G) = |G). But 
for any operator Y acting on to or less qubits we have 
(G|y|G) = 2~ m tr(y) because the density matrix of any 
to qubits is proportional to 1. Thus any stabilizer X £ S 
acts on at least to + 1 qubit. It means that |G) is an 
additive quantum code coding qubits into 2to qubits 
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with the minimal stabilizer weight m + 1 or greater. The 
results of the work |6) imply that for m < 15 such codes 
exist only for m = 1,3. 

The tensor powers of hexacode state |H® fe ) have the 
highest (to our knowledge) entanglement entropy for suf- 
ficiently small k. For example, the state |H® 4 ) has entan- 
glement entropy greater than determinant state |Det24,2), 
if one takes twenty-four qubits, see Table I. However we 
do not expect that |H) has the maximal entanglement 
entropy if all six qubit states could be considered. 



VII. CONCLUSION 

Although the generalization of entanglement entropy 
to multipartite case suggested in the present work looks 
rather natural, one faces a lot of difficulties while try- 
ing to compute entanglement entropy of some particular 
state. A progress can be achieved only if the state has 
some special properties or symmetry. 

For fixed n and d the maximal entanglement entropy 
S*(n, d) is rather close to the upper bound n\og 2 d, such 
that the ratio S*(n,d)/n\og 2 d approaches to one for 
fixed d and sufficiently large n. 

There is also one subtle point in the definition (^) con- 
cerning its stability under the extension of each party 
space of states. Suppose that each of the parties sharing 
the state |^) G adds new local degrees of free- 

dom to his part of the system thus extending his space 
to C D , D > d. The original space C d is somehow em- 
bedded to extended space C D and the original state \*ff) 
now is a vector from (C D )® n . There are two ways to 
compute entanglement entropy of \^): the parties may 
performe complete von Neumann measurements either in 
the original space C or in the extended space C D . Be- 
ing the functional of the state only, entanglement 
entropy should be the same in both cases. If S[^f] is in- 
deed invariant under such extensions, call |^) the stable 
state. In the bipartite case any state is stable, because 
S[$>] is invariant functional of single party density ma- 
trix, see (Q,|h. The determinant state |Det„) is also stable 
for any n. Indeed, the formula (^) which guarantees the 
equality S^Det^] = log 2 (n!) remains valid even if we al- 
low the states \<j>i) to be chosen from the extended space: 
the maximum is obviously achieved when all |</>i)'s belong 
to original space C n . However we can not prove that ar- 
bitrary multipartite state is stable, so this question is 
open. 

Also it is interesting to check whether generalized en- 
tanglement entropy is monotonic under local quantum 
operations for n > 2. 
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APPENDIX A: 

The purpose of this section is to prove the equal- 
ity (|2~4|). We will consider probability distributions p(x) 
of n classical bits, i.e. x — (x%, . . . , x n ), Xi = 0,1. 
By analogy with the set Pg we will consider the sets 
P^j of fc-uniform distributions of n bits. By definition, 
p(x) G P* iff any k of n bits have absolutely uniform 
distribution (i.e. if the sum of p(x) over any n — k bits 
is equal to 2~ fe ). For example, the set P™ consists of just 
one point — absolutely uniform distribution of n bits 
p(x) = 2~" while the set P° includes all possible n-bit 
distributions. If p{x) G Pj°j, consider its binary Fourier 
transform q(y): 



?(y) = E !ce B«(-i) a: ' 3/ p(^), 

p(x)=2- n ^ yeBn (-irvq(y). 



(Al) 



Here B™ denotes the set of all length n binary strings and 
x-y = J2"=i x iyi (mod 2). If y G B", denote wt(y) € [0, n] 
the number of T' in the binary string y. The definition of 
fc-uniformity can be rephrased in terms of Fourier com- 
ponents as: 



q(y) = if 1 < wt(y) < k, 

iff \ 

E, eB "(-irMy)>o. 



(A2) 



Note also that P^ is a convex set: if p'(x),p"(x) e P„ 
then for any a £ [0, 1] we have ap'(x) + (l— a)p"(x) G P*. 
It is known that the entropy is a concave functional, i.e. 

S[ap'(x) + (1 - a)p"(x)] > aS[p'(x)] + (1 - a)S\p"(x)}. 

It means that the minimum of 5 , [p(cc)] over all p(x) G Pj*j 
is achieved when p(x) is an extremal point of P„. Because 
P^ is specified by a finite number of linear equalities and 
inequalities, it has a finite number of extremal points. In 
principle we could find all them, compute S[p(x)\ for each 
point and choose the minimal value. However the set Pg, 
which we are interested in, has too many extremal points 
and such method doesn't work in practice. 

Instead we will proceed as follows. Consider six bit 
probability distribution p(x) defined by (p0|). An explicit 
verification shows that p(x) G Pg (as it should be, be- 
cause any measurement on the state |H) produces proba- 
bility distribution from Pg) and that 5[p(x)] = 4. It tells 
us that 
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inf S[p(x)} < 4. 

p(x)€P 3 6 



Also note that 



inf S\p(x)}< inf %(*)]. 



(A3) 



(A4) 



Indeed, take any p(x) G Pg and average out the sixth bit. 
Then, by definition, the distribution of the bits 1, 2, . . . , 5 
is 3-uniform: J2 Xe P( x ) ^ MS* ^ ne en tropy of the par- 
tial distribution can not exceed the entropy of the joint 
distribution, so that p(x)] < S[p(x)], which im- 

plies ( A4 ) . Now if we will manage to prove that 



inf S[p(x)}=4, 
p(z)ePj? 



(A5) 



then the work is done because (A3,A4,A5) imply (24) 



The remaining part of the text is the proof of equal- 
ity Q£§). 

The proof is based on extremal point analysis as it was 
suggested above. Let us parame terize p (x) G P| using 



its Fourier transform q(y), see (Al,A2). The nonzero 
components of q(y) are listed below: 



9(1,1,1,1,1) 
9(0,1,1,1,1) 

9(1,0,1,1,1) 



= 9, 9(1,1,0, 1,1) = 93, 
9i, 9(1, 1,1, 0,1) = 94, 
q 2 , g(l,l,l,l,0) = q 5 . 



(A6) 



Then p(x) can be written as: 
p(x) = (1/32){1 + (-l) wt (*> [q + J2 9,(-l) K! ]}• (A7) 

i=l 

The positivity constraint p(x) > 0, x G B 5 specifies the 
convex set P| in the space of q, qi, . . . , 35. If p(x) is an ex- 
tremal point of Pf then p(y) = for at least one y G B 5 . 
By the symmetry, we can assume that p(0, 0, 0, 0, 0) = 0. 
So to find all extremal points of Pf one suffices to find 
all extremal points of the convex set 



{p(x) G P5 : p(0,0,0,0,0) = 0}. 



(A8) 



In terms of variables q, 91, ... ,95 the set Pf is described 
by the following linear constraints: 



9 = -1 - Ei=i9i, 

ELi% > -1, 

9i + Qj < 0, 1 < i < j < 5, 



(A9) 



(fortunately, it appears that only part of inequalities 
p(x) > is independent). It is convenient to introduce 
one more auxiliary set Q defined as 

Q = {(9i,...,95) : 9*+9, <0, 1 < i < j < 5}. (A10) 

It is also a convex set. One can easily show that Q has 
only one extremal point 91 = • • • = 95 = and two 



types of one-dimensional edges with five edges of each 
type coming out from this extremal point: 



'i 

,(2) 



= {qi <0, 9j = if j 7^ i}, ie [1,5], 
= {9i > 0, qj = -qi, if j ^ i}, ie[l,5]. 



(All) 



The extremal points of P| are those extremal points 
of Q for which Ei=i Qi — ~1 anc l a l so the intersec- 
tions of one-dimensional edges of Q with the hyperplane 
Ei=i 9i = — 1- Summarizing, there are only eleven ex- 
tremal points of P| : 

1) 9 = -1, 9i = ••• = 95 = 0, 

2) g = 0, qi = -1, = if j 7^ i; i G [1,5], 

3) 9 = 0, q { = 1/3, q, = -1/3 if j ^ t; * G [1,5]. 

Here extremal points 2) and 3) represent the intersec- 
tions of and ef^ correspondingly with the hyper- 
plane Ei=i 1i = — 1 while 1) is the extremal point of 
Q. Substituting them into ( |A7| ) one can find the cor- 
responding distributions p(x). One can check that for 
extremal points 1) and 2) the probability p(x) takes only 
values and 1/16 thus having entropy iS[p(a;)] = 4. 
For extremal points 3) the probability takes the val- 
ues 0, 1/12, and 1/24. The entropy appears to be 
S[p(x)] = 17/6 + log 2 3 w 4.4. Thus for all extremal 
points of Pg we have iS[p(a;)] > 4 and for some extremal 
points S[p(a;)] = 4 which implies the equality (A5). 
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